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We apply the statefinder hierarchy plus the fractional growth parameter to explore the extended Ricci dark 
energy (ERDE) model, in which there are two independent coefficients a and j}. By adjusting them, we plot 
evolution trajectories of some typical parameters, including the Huhble expansion rate E, deceleration parameter 
q, the third- and fourth-order hierarchy Sj*' and 5^'* and fractional growth parameter 6 , respectively, as well as 
several combinations of them. For the case of variable a and constant /3, in the low-redshift region the evolution 
trajectories of E are in high degeneracy and that of q separate somewhat. Flowever, the ACDM model is 
confounded with ERDE in both of these cases. S 3 *' and S^'*, especially the former, perform much better. They 
can differentiate well only varieties of cases within ERDE except ACDM in the low-redshift region. For the 
high-redshift region, combinations jSj,*',e) can break the degeneracy. Both of jS 3 **, 6 ) and {S^’*,e) have the 
ability to discriminate ERDE with a = 1 from ACDM, of which the degeneracy cannot be broken by all the 
before-mentioned parameters. For the case of variableand constant a, and S^’*(z) can only discriminate 

ERDE from ACDM. Nothing but pairs jS^'^.e) and {S^*',e) can discriminate not only within ERDE but also 
ERDE from ACDM. Finally, we find that S 3 *' is surprisingly a better choice to discriminate within ERDE itself, 
and ERDE from ACDM as well, rather than S^*'. 


1. Introduction 

Data from a series of astronomical observations for more 
than a decade have shown that the universe is undergoing an 
epoch of accelerated expansion [1-4]. The most likely ex¬ 
planation for this cosmic acceleration is that the universe is 
currently being dominated by an exotic component, named 
dark energy (DE), which exerts repulsive gravity. To explain 
the origin and physical properties of dark energy, numerous 
theoretical/phenomenological models have been proposed [5]. 
Among the models, the most successful one is the ACDM 
model (which mainly includes the cosmological constant A 
and cold dark matter), because it is simple but could provide 
a very good fit to the observational data currently available. 
The cosmological constant is equivalent to the vacuum en¬ 
ergy density with w - -I. For a time-dependent equation- 
of-state parameter (EOS) w, there are lots of models, such 
as quintessence [6], Chaplygin gas [7], holographic dark en¬ 
ergy [8], and so on. 

In this paper, we study the model inspired by the holo¬ 
graphic principle of quantum gravity. The holographic princi¬ 
ple was enlightened by quantum properties of black hole [9, 
10] and later extended to string theory [11]. According to the 
work of Cohen et al. [12], when pde is taken as the quantum 
zero-point energy density caused by a short distance cut-off, 
the total energy in a region of size L should not be more than 
the mass of a black hole of the same size, i.e., L-'pde < 

The saturated form of this inequality, which is equivalent to 
the largest L allowed, leads to the energy density of the holo¬ 
graphic dark energy, pde = where c is a numerical 

constant introduced and Mp is the reduced Planck mass with 
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Mp = (SttG)”*. Eor the model setting, the choice of the in¬ 
frared (IR) cut-off L is very crucial. After the denial of the 
Hubble scale [13] and the particle horizon [10, 14] as IR cut¬ 
off for their failure to give rise to the cosmic acceleration, Li 
chose the future event horizon instead, getting the expected 
success [8]. But the adoption of the future event horizon in¬ 
dicates that the history of dark energy depends on the future 
evolution of the scale factor a(f), which violates causality [15]. 
Then the agegraphic dark energy model [15, 16] and the Ricci 
dark energy (RDE) model [17] emerged to avoid the viola¬ 
tion of causality. The former is characterized by the age of 
the universe as the length measure while the latter takes the 
average radius of Ricci scalar curvature as the IR cut¬ 

off. Further, the RDE model was extended to a more general 
form, liberating coefficients of the two terms, of the energy 
density [18, 19] 

Pde = 3MpV//"+m (1) 

where a and fi are constants to be determined and the dot de¬ 
notes a derivative with respect to time. For the Ricci-type 
holographic DE models are determined by a local concept 
of Ricci scalar curvature rather than a global one of future 
event horizon, they are naturally free of the causality problem. 
The model is called the extended Ricci dark energy (ERDE) 
model, with the special case a -2(3 the RDE model. 

With the increasing number of DE models, diagnostics aim¬ 
ing to differentiate them are needed. So far several methods 
have appeared. They are the well-known statefinder [20, 21], 
Om [22] and growth rate of perturbations [23-25]. The 
statefinder is a sensitive and robust geometrical diagnostic 
of DE, which uses both the second and the third derivatives 
of a{t). Recently, Arabsalmani and Sahni further extended 
the statefinder to higher-order derivatives of a{t), and called 
such a diagnostic “statefinder hierarchy” [26]. The statefinder 
diagnostic has been applied to various DE models [27-37], 
but sometimes we do need the diagnostic with higher-order 
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derivatives of a(t). For instance, when diagnosing the new 
agegraphic DE model, the original statefinder (second and 
third derivatives) cannot differentiate this model with different 
parameter values [36], but the hierarchy (further higher-order 
derivatives) is capable of breaking the degeneracy [38]. 

Here we study the ERDE model with statefinder hierarchy, 
supplemented by the growth rate of perturbations, to explore 
what the behaviors are like when ERDE takes different pa¬ 
rameter values and what the difference is between ERDE and 
ACDM. Necessarily referring to, in our previous work [35], 
we have diagnosed, with the original statefinders, the ERDE 
model both with interaction between DE and matter and not. 
The results therein seems satisfactory since there is no appear¬ 
ance of degeneracy for ERDE with various parameter values. 
But from the aspect of completeness of a theory, we neglected 
another evolution tendency of ERDE, mentioned in some pa¬ 
pers [39-41], that a of values larger than 1 enables the ERDE 
model to exhibit another orientation of evolution symmetrical 
to that plotted in Ref. [35]. As a matter of fact [40], a > 1 
makes ERDE behave like quintessence [ 6 ] (w > -1), while 
O' < 1 like quintom [42] (w evolves across the cosmological- 
constant boundary w = -1). We will expatiate on this theme 
in a later section. 

In Sect. 2, the ERDE model is exhibited. In Sect. 3, we 
introduce the diagnostic tools of statefinder hierarchy and 
growth rate of perturbations. Then ERDE will be explored 
in Sect. 4. Einally Sect. 5 gives the conclusion. 


2. The ERDE model 

We consider a flat universe with DE and matter, namely, 

3MIh^ ^ pa, + Pn,, (2) 

where pa, and pm are, respectively, energy densities of DE and 
matter, and pa, takes the form of ERDE described by Eq. (1). 
We get 


The fractional density and EOS of ERDE are given by 
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3. Statefinder hierarchy and growth rate of matter 
perturbations 

3.1 The statefinder hierarchy 


The primary aim of the statefinder hierarchy is to single the 
ACDM model out from evolving DE ones [26]. Eor it has a 
convenient property that all members of the statefinder hierar¬ 
chy can be expressed in terms of some elementary functions 
(like the deceleration parameter q, the EOS w or the fractional 
density Q), even the Chaplygin gas, interacting dark energy, 
and modified gravity models have already been explored in 
this way [43-45]. 

To review briefly, we just explain the primary principle of 
the statefinder hierarchy. Because in Ref. [26], the EOS of DE 
w in the hierarchy expressions is a constant, we here general¬ 
ize w to be time-dependent and have intensive interest only 
in the final expressions of statefinder hierarchy members in 
terms of elementary functions. Later we will see that these 
elementary functions are the fractional density and EOS of 
ERDE, already derived above, Eqs. ( 6 ) and (7). 

We Taylor-expand the scale factor a(t) around the present 
epoch to'. 


a(t) 
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where 


, i/2 ,, , j8dE^ 

(3) 

where E - H/Hq is the dimensionless Hubble expansion rate, 
X = \n a and the subscript “ 0 ” denotes present values of 
physical quantities. The solution of Eq. (3) is 
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n ~ 


a'") 


ne N; 


(9) 


a*'"'* is the nth derivative of a{t) with respect to time. The fa¬ 
miliar term A 2 = -q represents the deceleration parameter, 
while A 3 is the very original statefinder “r” [ 20 ]. A 4 was ever 
referred to as the snap “s” [46], while Aj the lerk [47]. Eor 
ACDM (w = -1), 




2a-3/3-2 




(4) 


where 


/o 


1 H- 


2a-3/3-2 


-a 


'm 0 > 


under the initial condition Eq — E(x - 0) - 1. 
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A 2 - I - 

A3 = 1, (10) 

32 

A 4 = 1 - etc, 

where Qm = |(1 -H g'), which means that for ACDM the ele¬ 
mentary functions are the deceleration parameter q or the frac¬ 
tional density parameter of matter Qm, because w is constant. 
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Then the statefinder hierarchy Sn can be defined as [26]: 

3 

S2 - A2+ 2 ^m, 

Si -Ai, ( 11 ) 

32 

S4-A4+—Q.m, etc. 

Comparing Eq. (11) with Eq. (10), one can get the essen¬ 
tial feature of this diagnostic that all the S„ parameters stays 
pegged at unity for ACDM during the entire course of cosmic 
expansion. 


^hIaCDM - 1- 


( 12 ) 



In fact, that is why are defined in this way, to differ from 
both other constant-w DE models and evolving ones. 

Remember that in Ref [20] there is a statefinder pair (r, s], 
where r is Sj and s = ^ belongs to the third 

derivative hierarchy and serves the aim of breaking some 
of the degeneracy present in r. To normalize letters of the 
alphabet, Arabsalmani and Sahni introduced a general pair 
[26] 


and 


5^') = A4 -h3 (1-H^), etc, 



- 1 


(14) 


where y is an arbitrary constant and the superscript “(1)” is 
used for discriminating not only between the original hier¬ 
archy S„ and but also between and its derivative 
Therefore, = [1,0] for ACDM and (r, s] is 

just ^ 3 ^^] with y = 3. In this paper we use only the 5$,*^ 
series as follows: 


FIG. 1: (color online). The evolution trajectories of the equation of 
state w versus redshift z of ERDE for variable a with fi = 0.5. Herein 
n„o = 0.27. 


where f{z) - dindjdina represents the growth rate of lin¬ 
earized density perturbations [25], 

f{z) - q;;(z), (19) 



3(1 - w)(l - Iw) 

+ --^( 1 -Qn,), 

125(l-fw) 


( 20 ) 


where w either is a constant, or varies slowly with time. 
Combine the fractional growth parameter e(z) with statefinder 
hierarchy to define a composite null diagnostic (CND): 
(5„, e] [26]. Eor ACDM, y 0.55 and e = 1 [25, 48], there¬ 
fore, (5„,e]-(l,l]. 


4. Exploring ERDE with statefinder hierarchy 
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where the prime denotes the derivative with respect to x = 
In a. 


3.2 The growth rate of perturbations 


The fractional growth parameter e(z) [23, 24] can supple¬ 
ment the statefinders as a null diagnostic as well, defined as 


e(z) := 


f(z) 


/acdm(z) 


(18) 


At first, it is necessary to clarify the different considera¬ 
tions between our previous work [35] mentioned above and 
this one, which can help to choose proper typical values of 
parameters here. In the former, we used an imposed condi¬ 
tion Wo = -1 to reduce by one the degrees of freedom for 
parameters a and /?, which were both born arbitrary. But in 
this paper, from the perspective of simulation in theory, we let 
w remain free so as to investigate the dependency level of the 
ERDE model upon a and fi by adjusting them. Therefore, we 
explore two cases. The first is adjusting a with a constant /?; 
the second is, whereas, adjusting /? with a constant a. 

For constant (3, we take (3 - 0.5 as a typical value, ap¬ 
proaching the observational constraints [17, 49, 50]. We fix 
f2m0 = 0.27 throughout the paper. Figure 1 exhibits the phe¬ 
nomenon we mentioned above, namely, “a > 1 makes ERDE 
behave like quintessence while a < 1 like quintom”. It is ob¬ 
vious that a plays a key role in the evolution of ERDE. When 
a > 1, the EOS evolves in the range of -1 < w < 0. When 
a < 1, the EOS evolves from the region of w > -1 to that 
of w < -1, i.e., the model exhibits a quintom-like evolution 
behavior. In particular, the boundary a - 1 can as well make 
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FIG. 2: (color online). The evolution trajectories of E, q, S^'' and 5^*', respectively, versus redshift z of ERDE for variable a with = 0.5, as 
well compared with the ACDM model. Herein = 0.27. 
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FIG. 3: (color online). The evolution trajectories of 5 3 ** and 5®, respectively, versus e of ERDE for variable a with fi = 0.5, as well compared 
with the ACDM model marked by a star. The dots represent the present values and the arrows indicate the directions of evolution. Herein 
Dmo = 0.27. 


the model behave like quintessence, but the universe will ulti¬ 
mately enter the de Sitter phase in the far future. It is neces¬ 
sary to emphasize that these features undoubtedly hold for the 
RDE model [39]. In our previous work [35], although (3 took 
0.3, 0.4, 0.5, and 0.6, under the condition of the present EOS 
value Wo = - 1 , all the values a obtained there were less than 
1. Just because of this, with no lack of universality, in this 
paper we explore the ERDE model comprehensively without 
missing any possibilities. Thus, around the boundary 1, we 
take a to be 0.8, 0.9, 1.0, 1.1, and 1.2. 


Eigure 2 shows the evolutions of the various-order deriva¬ 
tives of the scale factor a versus redshift z, from the hrst to the 
fourth, for the ERDE model. They are E, q, 5 3 *' and 5^ \ re¬ 
spectively, and they are meanwhile compared with the ACDM 
model. It can be seen that for E(z), according to Eq. (4), in the 
low-redshift region (z < 1 ) the curves of the model itself with 
various parameter values, even together with that of ACDM, 
are highly degenerate. Although the degeneracy is broken in 
the high-redshift region, however, the well-known observa¬ 
tional data are mainly within the low-redshift region, featur- 
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FIG. 4: (color online). A comparison of evolution trajectories of the 
fractional growth parameter e versus redshift z of ERDE for a = 1 
with p = 0.5 and that of the ACDM model. Herein 0^0 = 0.27. 


ing z < 1. For instance, for some supernova samples [51] 
the majority of the redshifts are in the range of z < 1 while 
only a few are in the range of a higher redshift, 1 < z < 1.4. 
Therefore, the current observations for Eiz) have not been of 
help so far. If the next generation Extremely Large Telescopes 
with high resolution would observe the high-redshift QSOs 
(2 < z < 5) [52], the evolution of E{z) may help effectively. 

For q{z), according to Eq. (15), in the low-redshift region 
the degeneracy that exists in the E{z) case is broken to some 
extent, but the trends of these curves are quite close to one 
another, including that for ACDM. Better exhibitions appear 
in the cases of and especially in the former. From 
the 53 ^^(z) plot in Fig. 2, the degeneracy is perfectly broken 
in the region of z < 1. The two symmetrical orientations of 
evolution due to different a, which has been concluded before 
from the w(z) plot of Fig. 1, appear apparently. When a > 1, 
namely w is always larger than - 1 , 5 3 *^ evolves decreasingly 
from 1 ; when a < 1 , namely w can evolve across - 1 , 5 3 *^ 
evolves increasingly from 1. But differing from 5 3 ^^ does 
not show the symmetrical aspect, but only in the same side, 
although the curves separate well in the low-redshift region. 
In both plots there is a common feature that when a - 1, 
getting 5 3 ^' = = 1, the same as that of ACDM. So in the 

and ^ 4 '(z) plots there are two shortcomings. One is a 
high degeneracy still existing in the high-redshift region. The 
other is the curves of a = 1 for ERDE superposing that of 
ACDM. In face of them, a single diagnostic of geometry fails 
to be achieved. Instead, we combine it with the fractional 
growth parameter, namely CND, trying to hnd a better way. 

Interestingly, when using CND, and of 

Eig. 3, we hnd that the degeneracy in the high-redshift re¬ 
gion can be broken clearly, especially ( 53 *', e] performing far 
better. To solve the second shortcoming, in both 5 3*^-6 and 
^ 4^-6 plots, ERDE with a = 1 exhibits a short line segment 
while ACDM is just a point [1,1]. We may see the reason at 
a glimpse of Fig. 4. That is to say, in the evolution history the 
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0.8 
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1.0 
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C (.ij 
^ 30 

2.088 

1.464 

1 

0.696 

0.552 

c(l) 

^40 

2.806 

1.492 

1 

1.043 

1.332 

eo 

1.0034 

0.9974 

0.9902 

0.9813 

0.9702 

30 



1.536 



^^40 



1.806 



Aeo 



0.0332 




TABLE I: The present values of statefinders and fractional growth 
parameter, and eo, and the differences of them, aSjq, 

and Aeo. For each case, aSjJ,’ = Sj^fmax) - Sjgfmin), AS^g = 
54 g(max) - 54 g(min) and Aeo = eo(max) - eo(min). 


a 



0.9 



p 

0.35 

0.40 

0.45 

0.50 

0.55 

o(l) 

^ 30 

1.580 

1.538 

1.499 

1.464 

1.433 

c(l) 

^40 

1.629 

1.578 

1.532 

1.492 

1.457 

eo 

0.9923 

0.9945 

0.9961 

0.9974 

0.9984 

30 



0.147 



40 



0.173 



Aeo 



0.0061 




TABLE II: The present values of statefinders and fractional growth 
parameter, Sjg, and eo, and the differences of them, aS ^J,*, AS^jg 
and Aeo. For each case, ASjg = Sjgfmax) - Sjgjmin), AS^g = 
540 (max) - 54 Q(min) and Aeo = £o(max) - eo(min). 


fractional growth parameter e(z) becomes closer and closer to 
1 from past to present, but not equal to yet. Since the present 
values of physical parameters are significant in the research 
of cosmology. Table I shows the present values of parameters 
‘^ 30 ’ ‘^40 ^ 0 , and the differences of them, for each case 

aS 3 g = S 3 Q (max) - S 3 ^ (min) and the same way for AS 4 Q and 
Aeo. We can see a^^q > A^^q, namely the fourth derivative 
of a, compared with the third derivative, alleviates the degen¬ 
eracy of present values. But even so, the comparison of either 
the S 3 *^-z and S 4 '-z plots in Fig. 2 or the S 3 '^-e and S 4 *-e plots 
in Fig. 3 shows S 3 ^^ to be performing much better during the 
evolution process than 54 *. This indeed violates our habitual 
judgement: that the higher the order of derivative is, the better 
the diagnostic performs [38, 43-45]. 

For constant a, although it can be either larger or smaller 
than 1 , which leads to reverse orientations of the evolution, 
we only take a = 0.9 as a typical value in this exploration, 
according to the best-fit values for a from some of the recent 
constraints [49, 50, 53]. To obtain feasible evolutions, we take 
0.35, 0.4, 0.45, 0.5, and 0.55 for (5. Likewise, we explore the 
four parameters E, q, 5 3 *^ and ^ 4 ' in Fig. 5 first, as well as 
make a comparison with the ACDM model. We can see that 
the ERDE model is insensitive to parameter (3. Eor the hrst- 
and second-order hierarchy Eiz) and q{z), high degeneracy 
appears, even together with ACDM. In the third- and fourth- 
order cases, the ERDE model itself with various parameter 
values highly degenerates, but the ACDM model can be dis¬ 
criminated perfectly from ERDE in the low-redshift region. 
Then let us observe the CND of Fig. 6 . The and ^ 4 '-e 
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FIG. 5: (color online). The evolution trajectories of E, q, and respectively, versus redshift z of ERDE for variable j 8 with a = 0.9, 
as well compared with the ACDM model marked by a star. The dots represent the present values and the arrows indicate the directions of 
evolution. Herein Gnjo = 0.27. 
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FIG. 6 : (color online). The evolution trajectories of 5 3 ** and 5respectively, versus e of ERDE for variable p with a = 0.9, as well compared 
with the ACDM model marked by a star. The dots represent the present values and the arrows indicate the directions of evolution. Herein 
= 0.27. 


curves look as good as the above case of j3 - 0.5 in Fig. 3. 
In both plots the evolution trajectories separate quite well, but 
the combination of e) is slightly better than e) be¬ 
cause of the more legible separation in-between curves in the 
high-redshift region. In the same way we show in Table II the 
present values of and eo, and the differences of them 

for a - 0.9. Likewise the relation AS > AS j'q demonstrates 
once again that the fourth-order hierarchy can help to allevi¬ 
ate the degeneracy of present values when compared with the 


third one. But for the same reason as of the comparison of 
S^^^-e and 5^'^-e plots in Fig. 6 , we hnd for the ERDE model 
e) is a more efficient parameter pair of diagnostic than 
{54 \e), which is already concluded in the above-mentioned 
case of p - 0.5. 
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5. Conclusion 

In this paper, we explore the extended Ricci dark energy 
model with statefinder hierarchy supplemented by the growth 
rate of perturbations. Since in ERDE there are two indepen¬ 
dent variables a and J3, we just adjust them, respectively, leav¬ 
ing other parameters fixed, for the sake of investigating the 
effects of a and jS on this model. Eirst, a feature of the holo¬ 
graphic Ricci-type dark energy models is corroborated again, 
namely, a > 1 makes them behave like quintessence while 
a < 1 like quintom. For the ERDE model with fi — 0.5, letting 
a vary around 1, we conclude that the evolutions of the Hubble 
expansion rate E are in high degeneracy in the low-redshift re¬ 
gion of z < 1 ; but because the observational data come mainly 
from the low-redshift region z < 1 , the broken degeneracy in 
the high-redshift makes no sense. The evolutions of deceler¬ 
ation parameter q do degenerate no more in the low-redshift 
region of z < 1. However, for both E and q plots, the evo¬ 
lution of ACDM cannot be singled out from in-between with 
great ease. The situations of 5 3 *^ and which contain the 
third and fourth derivatives of the scale factor, respectively, 
turn out to be better, evolves with respect to redshift z 

along two orientations symmetrical to each other on the basis 
of different a in the region of z < 1. When a > 1, it evolves 
decreasingly from 1 ; when a < 1 , it evolves increasingly from 
1 . 54 ' although seems featureless by contrast with by 

a comparison of the present-value differences of the param¬ 
eters 5 3 *' and = 1.536, = 1.806), we see 

that 5^^ is capable of alleviating the degeneracy existing in 
other statefinder parameters for ERDE. There are also two un¬ 
solved problems that high degeneracy still exists in the high- 
redshift region, and with a = 1 are degenerate with ACDM. 
As for them, the combination of statefinder hierarchy S „ and 
fractional growth parameter e (CND) can help. In both 
and S plots, the degeneracy in the high-redshift region is 


pretty broken and ERDE with a = 1 exhibits a short line seg¬ 
ment, but ACDM exhibits just a point {1,1). Nevertheless, the 
fact that the S 3 '^-z and S planes feature a more legible and 
regular sight of evolution due to a than that of the S ^ ^-z and 
S 4 planes, reveals that the third-order statefinder hierarchy 
5 3 '^ makes more sense for ERDE than 5^ ^ does. 

For the ERDE model with a = 0.9, we find by contrast 
with a, p has a weak influence. The evolution trajectories of 
E, q, 53 *' and 5^' with respect to redshift z are in high de¬ 
generacy within the ERDE model. The degeneracy of ERDE 
with ACDM still exists for E(z) and ^(z), but it is perfectly 
broken in the low-redshift region for both ^ 3 *^^) 

As for the high-redshift region, the use of CND can break the 
degeneracy there, especially the { 53 *\ e) pair performs more 
efficiently than the pair, although A^jg = 0.147 < 

A^^g^ = 0.173. ACDM can also be discriminated from ERDE 
by CND. As a consequence of all the materials studied above, 
we find that, although the higher-order statefinder hierarchy, 
even with the growth rate of perturbations, can differentiate 
the ERDE model itself with various parameter values and also 
from the ACDM model, there is the interesting discovery that 
the third-order hierarchy of statefinder is really a better choice 
than the fourth-order hierarchy for the ERDE model. 
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